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Abstract — We study the question of determining locations of 
base stations that may belong to the same or to competing 
service providers. We talie into account the impact of these 
decisions on the behavior of intelligent mobile terminals who 
can connect to the base station that offers the best utility. The 
signal to interference and noise ratio is used as the quantity that 
determines the association. We first study the SINR association- 
game: we determine the cells corresponding to each base stations, 
i.e., the locations at which mobile terminals prefer to connect to 
a given base station than to others. We make some surprising 
observations: (i) displacing a base station a little in one direction 
may result in a displacement of the boundary of the correspond- 
ing cell to the opposite direction; (ii) A cell corresponding to a BS 
may be the union of disconnected sub-cells. We then study the 
hierarchical equilibrium in the combined BS location and mobile 
association problem: we determine where to locate the BSs so as 
to maximize the revenues obtained at the induced SINR mobile 
association game. We consider the cases of single frequency band 
and two frequency bands of operation. Finally, we also consider 
hierarchical equilibria in two frequency systems with successive 
interference cancellation. 



I. Introduction 

In this paper we study some hierarchical decision making 
problems arising in the uplinks of cellular networks. We 
first address the problem of association: given multiple base 
stations (BS) capable of providing services to a mobile located 
at a given point in the region of operation, to which BS should 
the mobile connect? This is studied in a non-cooperative 
context where each mobile connects to the BS that provides it 
with the best signal to interference and noise ratio (SINR). The 
associations determine the cells corresponding to each BS. We 
characterize the nature of cells as a function of BS locations. 

We then consider the problem of determining the locations 
of base stations, taking into account the behavior of the 
mobiles that will be induced by the location decisions. We 
study cases where the BSs cooperate (e.g., they belong to 
the same service provider) and those where they compete 
with each other The latter scenario results in a location game 
between the BSs. 
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Related work: Plastria H] presented an overview of 
research on locating one or more new facilities in an envi- 
ronment where competing facilities already exist. Gabszewicz 
and Thisse ||2| provided another general overview on location 
games. Mazalov and Sakaguchi [S) and references therein 
studied competition over prices of goods between facilities 
that have fixed positions. They then derived the equilibrium 
allocation of customers. Such games, as well as hierarchical 
games in which firms compete for location or over prices 
which then determine the customer-allocation equilibrium, 
were introduced by Hotelling ||4] in 1929. When considering 
such games over a finite line segment with two firms, the 
models under appropriate conditions give rise to a partition 
of the segment into two convex subsegments or "cells" as 
introduced in our context. 

An interesting difference between the settings above and our 
setting, which is also defined on a finite line segment, is that 
in our case more complex cells are obtained at equilibrium. 
This is due to the difference in the cost structure in the cellular 
context. Hotelling H considered a general cost related to the 
distance between the customer and the firm it chooses; this 
cost however depended only on the distance and not on the 
actual location of the firm. This does not hold in our case: 
the throughput of the mobile depends on the interference at 
the base station which in turn depends on the location of the 
base station. We finally note that in our model, the power of a 
mobile, which can be considered as the "cost", is fixed, while 
it attempts to maximize the utility, i.e., throughput. 

Aram et al. ||5| study coalition based joint resource pro- 
curement and resource allocation in wireless networks using 
the framework of cooperative game theory. They consider a 
set of operators and customers with predetermined customer- 
operator associations. All the operators together place base 
stations, procure spectrum, and allocate channels to the com- 
mon pool of customers. Doing so is shown to be optimal 
even when the operators are selfish. However, it is assumed 
that operators can divide the aggregate earned utility in any 
arbitrary way. In another work, Aram et al. [0 extend the 
analysis to nontransferable utilities. 

After the appearance of our initial work ||2], Ramanath 
et al. |[8l have studied the joint placement of two base stations 
that use same frequency. The utilities of BSs come from the 
notion of a-fairness. They also consider a multicell scenario 
where each cell has one BS, and all BSs operate on the 
same frequency. The goal is to optimally place the BSs in 
their respective cells so as to maximize their a-fair utilities. 
Silva et al. ||9l consider an association problem where mobiles 
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have hard average throughput constraints, and the objective 
of optimal association is to minimize the aggregate power 
consumption of all the mobiles in the network. They also study 
the downlink scenario under the assumption that neighboring 
BSs operate in orthogonal channels. Again the objective is to 
find an optimal association that minimizes the aggregate power 
consumption of the BSs. Kasbekar et al. flOl consider a joint 
problem of mobiles' association and charging and spectrum 
leasing by service providers. 

A description of the model studied in this paper and the 
notation used can be found in Section and Appendix [A] 

Our contributions: First, we consider a code division 
multiple access (CDMA) system where BSs perform single 
mobile decoding. We derive analytical expressions for the 
cell boundaries in the case where BSs are on the same 
frequency band (Section IIII-Al i. This allows us to study the 
geometric properties of cells as a function of the locations 
of the BSs. We then study the hierarchical equilibrium in the 
combined BS location and mobile association problem, i.e., 
we determine where to locate the BSs so as to maximize 
the revenues obtained on the induced SINR-based mobile 
association game (Section IIII-BI ). We also do the analogous 
analyses for the case where BSs are on different frequency 
bands (Section IIVI ). Subsequently, we consider BSs capable 
of successive interference cancellation (SIC) decoding. After a 
discussion on the association problem and the single frequency 
band case (Section 0, we analyze the case of different 
frequency bands and give a complete characterization of the 
resulting equilibria (Section IVII i. 

While the main body of our work assumes mobiles placed 
over the line-segment [—L, L] C M, we discuss the extension 
to two dimensional deployments in Appendix iHl 

II. The model and notation 

Our focus is on communication in the uplink direction, 
i.e., from mobiles to BSs. Mobiles and BSs lie in the two- 
dimensional space R^. A large number of mobiles are placed 
uniformly over the segment [—L, L] on the first of the coordi- 
nate axes. The fluid approximation is obtained for the infinitely 
large population of mobiles. For details see Appendix lAl There 
are two BSs, BS 1 and BS 2, located at (a-i, 1) and (x2, 1), 
respectively (say on the top of a flat building whose height is 
one unit). BSs cooperate if they belong to the same operator, 
and compete if they belong to different operators. We allow 
placements of BSs outside the area where mobiles exist, i.e., 
Xj < —L and Xj > L are allowed for j = 1,2. In the 
following, we use only the first coordinates to specify locations 
(with the understanding that second coordinates are in case 
of mobiles, and 1 in case of BSs). 

Transmitters are point sources radiating in two-dimensional 
space with circular wavefronts (respectively, three-dimensional 
with spherical wavefronts). We consider a power law path loss 
model with exponent a; i.e., the power from a radio transmitter 
attenuates as distance raised to the power a (see Appendix lAl). 
A mobile located at y has a channel "gain" of [{y — xj)"^ + 
^~a/2 -Qg J ^jj jjjobiles are assumed to transmit at a power 
such that the power density along the line is unit power per 



unit length. Thus the total transmitted power is 2L. Thermal 
noise at the BSs is assumed to be Gaussian with noise variance 
cr^ per sample. 

At any time, each mobile is associated with exactly one 
BS. Let Aj C [— L, L] be such that the mobiles in Aj are 
associated with BS j. Aj will be called cell j. The utility of 
a mobile at y is assumed to be a nondecreasing function of 
the SINK density at y, as seen at the BS to which the mobile 
is associated. The SINR density depends on the interference 
model under consideration, which we discuss next. 

A. Interference models 

Mobiles that connect to a particular BS may or may not 
cause interference to the other BS depending on whether the 
BSs operate on the same or different radio frequency (RF) 
bands. We consider both the cases in this paper. The case in 
which the same frequency band (channel) is used at both the 
BSs occurs if the wireless network operates in an unlicensed 
band; in such a case, BSs belonging to different networks (or 
providers) may use the same RF band. We call this the single- 
frequency case. If the wireless network operates in licensed 
RF bands, two neighboring BSs would operate in disjoint RF 
bands. We call this the two-frequencies case. We now discuss 
the useful power collected and interference seen at a BS in 
the single- and two-frequencies cases. For this purpose, it is 
useful to define the following functions. Define 

g(y):=[l + yT"^\ (1) 

For a set C [—L, L] and candidate BS location x, define 

E{x,S) I 9{y~x) dy. (2) 
Js 

and E°{x) := E{x, [—L,L]). The dependence of g, E, and 
E° on a is understood. 

In the following we consider a CDMA system where BSs 
perform single mobile decoding, i.e., while decoding any 
mobile's signal, they treat all other mobiles' received signals as 
interference. Subsequent descriptions of SINR-equilibrium and 
hierarchical equilibrium are also for such a system. Analogous 
notions for SIC decoding are defined in Sections |V] (single- 
frequency case) and |VT] (two-frequencies case). 

1) The single-frequency case: In this case, power from 
all the mobiles is received at both the BSs. The total re- 
ceived power at BS j located at Xj is therefore given by 
E{xj,[—L,L]) ~ E°{xj). All of this received power will 
clearly be interference to a mobile at y because the mobile's 
own contribution to this is infinitesimal. 

With this interference interpretation for E°{xj), we now 
highlight some of its properties. It is straightforward to see 
via change of variables that 

/L — x />arctan(L — a;) 

g{y)dy= / (cose)"-^ (3) 

-L — x J a.rctSin{ — L—x) 

Closed form expressions are available for E° when a takes 
integer values. In particular, for a = 2 we get 

E°{x) — arctan(L — x) + arctan(_L + x), (4) 
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and for a = 1 we get 



E°{x) — arcsinh(L ~ x) + arcsinh(i + : 



(5) 



The above expressions motivate the following definition of the 
a-parametric function 



arctanc (x) 



g{y) dy, x G 



Then clearly arctaiiQ. (•) is an odd functiorQ that is increasing, 
differentiable with derivative g, and sigmoidafl We may 
therefore write the received power at location Xj (and therefore 
the interference in the single-frequency case) as 



E°{x.j) = arctaiiQ, [L — Xj) + arctaiic {L 



(6) 



The following is a useful property of E°. 

Proposition 2.1: E° is an even function with a unique 
maximum at 0. Moreover, i?°(|a;|) monotonically decreases 
with |.t|. 

Proof: See Appendix iB] ■ 
2) The two-frequencies case: Unlike the previous setting, 
in which the two BSs operate on the same RF band, in 
the two-frequencies case the total interference at each BS 
depends on the association decisions of the mobiles. Indeed, 
the interference power at BS j is the total power received at 
that BS from all mobiles that actually associate with it. The 
total received power at BS j is thus given by E{xj,Aj). 

For example, suppose Ai := [—L,0\ and A2 := {d,L\ 
denote the two cells for some 9 G [— L, L]. Then the inter- 
ference power at BS 1 is E{xi^Ai) = arctana (6* — xi) — 
arctaiiQ, {—L — xi). The expression for i?(.T2, A2) is obtained 
analogously. 

B. SINR-equilibrium association 

We shall first consider the case in which the BSs' locations 
are fixed, and each mobile has the option of associating with 
one of the BSs. The continuum of mobiles constitute the 
players in this association game. 

Consider a mobile at location y. Its utility is a nondecreas- 
ing function of the throughput density at y (see Appendix 
IaI i. The throughput density at y increases linearly with SINR 
density. Thus, this mobile chooses a BS that yields the higher 
SINR density at y. Let Ij be the set of interferers as seen at 
BS j. If a mobile at point y is associated with BS j, the SINR 
density for this mobile is 



SINR(2/,x,,/,) 



(7) 



E{x,J,) + <7^- 

A mobile at y e [—L, L] will therefore prefer to associate with 
BS 1 if SINR(y,a;i,/i) > SmR{y,X2, h)- 

We observe that in the single-frequency case, Ij = [—L, L]. 
Thus, the SINR density for a location, as seen at BS j is 
fixed. However, in the two-frequencies case, Ij ~ Aj,j = 1,2. 

'it is an odd function because arctaria (—a;) = — arctan^ (x). 

function is sigmoidal, if it is non-decreasing, concave to the liglit of a 
particular point called the inflection point and convex to its left. The second 
derivative of arctan^ (x) = g'{x) = — + x^]"'^"*""/^). The inflection 
point for arctan^ (■) is therefore 0. 



Hence, the SINR density for a location, as seen at BS j is a 
function of the cell Aj. 

Definition 2.1 : The cell partition {Ai,A2) is said to 
be an SINR-equilibrium if the following holds: y £ 
Ai if SINR(?/,a;i,/i) > SINR(y, 0:2, 12) and only if 
SINR(y,xi,/i) > SINR(2/,a;2,/2). If SmR{y,xi,h) = 
SINR(2/,a;2,/2), y e Ai or A2 arbitrarily. 

Remark 2.1: This definition of equilibrium is similar to 
the Wardrop equilibrium in road traffic [ITTI . or the Nash 
equilibrium in population games lfT2l . Note, however, that in 
Wardrop equilibrium the utility of choosing a resource (a BS 
in the present problem) depends on the set of users that make 
the same choice through their total "number" (their fraction 
or their mass). Extensions of the Wardrop concept exist to the 
case where there is a finite number of user classes and the 
utility of using a resource for a user in a given class depends 
on the amount of users of each one of the classes who use that 
resource Iil3i . In our problem, however, there is a continuum 
of classes corresponding to the locations of the mobiles. 

C. Hierarchical equilibrium 

We shall also consider placement of BSs taking into account 
the SINR-equilibrium that follows when mobiles associate to 
maximize their SINR density. The two BSs play a location 
game: BS j decides to place itself at {xj , 1) where Xj £R,j = 
1, 2. The utility of a BS is a monotone function of the aggregate 
throughput of all the mobiles associated with it. Since the 
throughput density at location y increases linearly with SINR 
density, we may simply set the integral of SINR density over 
the cell of a BS as its utility. Thus for BS j with cell Aj and 
interferers Ij, the utility is 



SINR(y,x,,/,) dy^]- 



gjy-xj) dy 
E{xjjj) + a'^ 



Once the BSs choose their locations, Aj, Ij, and thus the utility 
of BS i are determined by the association game played by the 
mobiles. We thus have a Stackelberg-like game llT4l with the 
lead players being the two BSs (who may either cooperate or 
compete) and the followers the continuum of mobiles (who 
compete to maximize their respective SINR densities). We 
refer to this as the hierarchical equilibrium problem. 

III. CDMA: The Single-Frequency Case 

A. SINR-equilibrium association 

We begin by providing closed form expressions for cell 
boundaries in the SINR-equilibrium (see Definition 12.1b . De- 
fine the ath root of the ratio of the net interferences (including 
thermal noise) at the two BSs to be 



Ba{xi,X2) 



E°{xi)+a^ 

E°{X2)+(J^ 



l/a 



We start by considering symmetric placements of BSs: |xi| = 
1x21 which implies Ba{xi,X2) = 1. If xi = X2 BSs are 
indifferent to all the mobiles from the point of view of 
SINR density. Hence {A, [-L, L]\A) for all A C [-L, L] 
are SINR-equilibrium association profiles. If xi = —X2 7^ 0, 
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mobiles associate with the BS which is closer. Hence either 
([-L,0],(0,L]) (if xi < 0) or ([0, L]), [-L, 0)) (if xi > 0) 
is the unique SINR-equilibrium. To study the asymmetric 
scenarios we assume, without loss of generality, that BS 2 
is located closer to the origin than BS 1, i.e., \xi \ > \x2\ > 0. 
On account of Proposition 12.11 we have Ba{xi,X2) < 1. 

Proposition 3.1: Let BS 1 be located at xi and BS 2 at 
X2 where |a;i| > \x2\ > 0. The set of mobile locations that 
connect to BS 2 is nonempty only if 



a = 2 



a = 1 



\X2 - Xi\ 



Ba{xi,X2) 
l-Bl{xi,X2) 



> 1. 



(8) 



If the inequality holds strictly then the set of locations that 
connect to BS 2 is given by the interval^ 



X2 - XiBl{xi,X2) 
1 - BI{X1,X2) 



+ (^-Vt^ - 1, - 1^ 



Proof: Mobiles that have a higher SINR density at X2 
will connect to BS 2, i.e., y E A2 if 

[(y-X2)2 + l]-"/2 ^ [{y - f + l]--/^ 



E°{X2) + 0-2 

which is equivalent to 

{y-x2f + l< {{y 



E^ixi) + 



+ Bl{xi,X2). 



As B'^{xi,X2) < 1, the above inequality holds when a 
convex quadratic function of y is strictly negative. The positive 
discriminant condition straightforwardly yields that the set 
connecting to BS 2 is nonempty only if (O holds. The roots 
of the convex quadratic equation are given by the ends of 
the specified interval. Since the convex quadratic function 
is strictly negative in the interval between the roots, all the 
mobiles in this interval have higher SINR densities at BS 2. 

■ 

When 1x2 1 > |xi| > 0, the roles of BS 1 and BS 2 are 
switched: BS 1 sees more interference, its cell Ai may be 
empty, and when nonempty, Ai is an interval. 

We provide numerical results to illustrate some surprising 
features of the SINR-equilibrium that distinguish this from 
other association games (e.g., ID, ifTSl ). We set L ~ 10 (so 
that mobiles are concentrated over the interval [—10, 10]) and 
the noise parameter a = 0.3. We place BS 1 at one of the 
fixed locations xi where xi = —10,-5,-2,0. For each of 
these, we vary the location of BS 2 from X2 = to X2 = 
30 (see Figure The left column of plots corresponds to a 
path loss exponent a ~ 2 and the right one to a = 1. The 
equilibrium sets Aj turn out to have the form Ai = [01,62], 
A2 = [-L, 0i) U (02, L] for xi = 0, -2, and Ai = [-L, 62], 
A2 = {62,14 for xi = —5, —10. The top (respectively bottom) 
row of plots depict the threshold 62 (respectively 6*1) as a 
function of X2- See the following for more details. 

1) Observations: 

^The notation a + {b, c) is short for the interval (a + b, a + c). 
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Fig. 1. Single-frequency case; SINR-equilibrium: Thresholds determining 
the cell boundaries as a function of the location of BS 2 for various locations 
of BS 1. The path loss exponent is 2 in the figures on the left and 1 in those 
on the right. 



a) Non-convex cells: For all the locations of BS 1, xi = 
— 10, —5, —2, 0, mobiles in (^2, i] have a better SINR density 
at BS 2. Let us concentrate on the curves corresponding to 
xi = —2 in Figure [T| When BS 2 is located sufficiently far 
to the right of the origin, the interference at BS 1 is large 
compared to that at BS 2 (see Proposition 12. It . Thus, mobiles 
sufficiently far away and to the left of BS 1 (those in [—L, 9i)) 
also have a better SINR density at BS 2 despite BS 2 being 
the farther BS. Thus, in this case, A2 = [-L, di) U (6*2, L], a 
non-convex set. A2 is similarly non-convex when xi = and 
X2 is sufficiently far to the right (or left). 
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Fig. 2. Single-frequency case; SINR-equilibrium: Upper cell boundary of 
cell A\ as a function of the location of BS 2 for various locations of BS 1. 

h) Non-monotonicity of the cell boundaries: We observe 
a surprising non-monotonicity of the threshold 62 as a function 
of the location X2 of BS 2. 6*2 first increases with X2 until about 
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X2 ~ 8, then it decreases with X2 until around X2 = 14; finally, 
for larger X2, O2 again increases. Analogous non-monotonicity 
is observed in 61 too (in the the curves corresponding to xi — 
-2,0 and a = 2). 

The dashed line in Figure [2] shows a zoomed-in view of the 
xi = —10 case of the top-left plot of Figure [T] The threshold 
62 increases beyond until X2 is about 8 units to the right of 
the origin, and then returns to when X2 = 10- This can be 
understood as follows. Clearly, for X2 = 10 the interferences 
at both the BSs are the same; hence 02 = 0, the midpoint. Now 
imagine moving BS 2 a little to the left (i.e., decreasing X2)- 
Now \x2\ < \xi\. Thus, from Proposition 12. 11 the interference 
E°{x2) at BS 2 is larger than E°{xi), the interference at BS 1. 
This makes it advantageous for mobiles a little to the right of 
the origin also to associate with BS 1; hence 62 increases as 
X2 decreases from X2 = 10. Further decrease in X2 brings 
BS 2 closer to mobiles on the negative .x-axis, thus ultimately 
causing 82 to return to 0, and even cross below 0, as X2 
decreases further As X2 increases beyond 10, the interference 
perceived by it decreases, thus making it advantageous for 
mobiles a little to the left of the origin also to associate with 
BS 2; hence 62 decreases as X2 increases beyond X2 ~ 10. 
Once BS 2 is moved far from the region where the mobiles 
exist, the signal power to X2 becomes smaller, and association 
with BS 1 becomes increasingly better for mobiles to the right 
of the origin, causing 6*2 to increase. 

The top row of plots in Figure [T] suggests that 62 is perhaps 
monotone in the position of BS 1. But this is not true because 
a closer look at the 6*2 curves in Figure |2]for xi = -10,-8 
shows that they cross each other several times. 

2) Discussion: The form of equilibria displayed in the 
SINR-association examples is unusual in the class of location 
games. The reason for the unusual features lies in the SINR 
criterion, as we describe now. 

• If a mobile is very close to a BS, path gain from the 
mobile to the BS will be very high. Thus, the mobile 
connects to this BS, even if the interference suffered by 
this BS is relatively higher 

• If a mobile is located sufficiently far from both BSs, then 
the relative difference in the powers received at the BSs 
will be small. Thus the mobile will prefer to connect to 
BS that suffers from less interference. 

• If a mobile is at moderate distance from both the BSs, 
it takes into account both the factors (i) path gains to 
the BSs and (ii) interferences suffered by the BSs, while 
making its association decision. 

B. Hierarchical equilibrium 

Single base station: Suppose there is only one BS. Given 
that the interference is maximum at the origin and decreases 
monotonically with distance from the origin, where should it 
be placed to maximize utility? The utility of the BS, when 
placed at x, is given by 

1 gjy - , ^ 1 E-{x) 

2 J_L E''{x)+a^ ^ 2E"{x)+a^ 



which is maximized when E°{x) is maximized, i.e., at x = 0. 
Despite the high interference, the origin is the best location to 
maximize the utility given the nature of the utility function. 



0.5 




0.1 



-^0 -10 10 20 

^2 

Fig. 3. Single-frequency case; Utility of BS 2 as a function of its location 
when we position BS 1 at xi = —10. 

Two base stations, utility behavior: Figure |3] portrays 
the effect of cell boundaries on utility. BS 1 is located at 
xi = —10. The utility of BS 2 as a function of its location 
is then plotted. For X2 <C —10, mobiles at the farthest right 
connect to BS 2. For X2 approximately in the interval -11 to 
-10, an interval in the left most extreme of [—10, 10] also 
joins BS 2 so that BS 2 cell partition is a union of two 
intervals. Consequently, the utility increases in this interval. 
At 2:2 = —10 a sudden transition occurs where all nodes in 
the middle interval switch to BS2, and hence the discontinuity 
in the utility. For — 10 < a;2 < 10 nodes in an interval join 
BS 2, and this eventually becomes a half infinite line with 
boundary moving to the right as X2 becomes large. 

1) Two cooperating base stations: We now consider opti- 
mal joint placement of two BSs to maximize the sum utility. 

It can be shown that in a hierarchical optimal configuration 
the two BSs are placed on the opposite sides of the origin. 
Furthermore, the corresponding cells are of the form [—L,a] 
and (a, L], and are characterized by a single parameter a. See 
Appendix ICl for justifications. While the exact characterization 
of a remains open, simulations indicate that sum utility is 
maximized when a = and —xi = X2, i.e., the BSs 
are equidistant from the origin. We call such a placement 
as symmetric. The SINR-equilibrium cells under symmetric 
placement are [—L,0] and {0,L]. 

Remark 3.1: The optimal configuration should also be an 
SINR-equilibrium association profile implying that for a user 
at a, SINR densities seen at BS 1 and BS 2 must match, i.e., 

gja-xi) ^ g{a-X2) 
E°{xi) + a'^ E°{x2) + (t'^' 

Note that the symmetric configuration with boundary a = 
and —xi = X2 satisfies this condition. 

Figure |4] depicts the utility obtained by each BS for sym- 
metric placement —xi = X2 = a:, as a function of x. We see 
that the origin and the extreme points (at distance 10 from 
the origin) are suboptimal locations. We also observe that the 
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Fig. 4. Single-frequency case; cooperating BSs; symmetric placement of 
BSs at ±x: Plot of utility obtained by each BS vs. x; here L = 10, a = 2 
and a takes the values 0.4, 1, and 2. 



Fig. 6. Single-frequency case; non-cooperating BSs: Utility of BS 2 as 
a function of its location when we position BS 1 at xi where xi = 
-2,-5,-8,-10. 



performance close to the optimal location is quite robust to 
perturbations of BS locations. 

lOi ' ' ' 1 



9 




^0 0.5 1 1.5 2 
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Fig. 5. Single-frequency case; cooperating BSs; symmetric placement of 
BSs: Optimal symmetric distance from origin for two BSs, as a function of 
thermal noise standard deviation. 

Further experimentation revealed that, as a is increased, 
optimal distance of the BSs from the origin decreases (see 
Figure |5]l. As rr oo, the optimal symmetric locations of the 
BSs converge to -5 and 5. This is expected because at very 
large a, interference does not play any role, and the BSs should 
be placed to maximize the total power collected from the 
respective cells, E{x, (0, L])+E{—x, [— L, 0]). Proposition l2.1l 
says that this is maximized by choosing x and —x to be the 
mid-points of the respective intervals, i.e., x = L/2, which is 
5 in our example. 

2) Two non-cooperating base stations: We now consider 
a non-cooperative game between the two BSs. The BSs act 
simultaneously and pick their locations to maximize their 
respective utilities. 

Figure |6] has on the horizontal axis the location of BS 2 and 
on the vertical axis the utility it achieves. The figure is obtained 
for L = 10, a = 0.3, a ~ 2. There are four curves that 
correspond to four locations of BS 1: xi ~ —2, —5, —8, —10. 
From these curves, one can conclude that the utility of BS 2 



is quite robust to placement errors around the best response 
location, for the indicated values of BS 1 locations. 

20| ' ' ' 71 
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0^ • • • 1 

5 10 15 20 

-x^: Distance of BS1 from origin 

Fig. 7. Single-frequency case; non-cooperating BSs: The best response of 
BS 2 when BS 1 is at a distance indicated by abscissa to the left of the origin. 

Figure I2] shows the best response of BS 2 to a BS 1 location. 
BS 1 is moved along the segment to left of the origin. In 
the figure the horizontal axis is — .ti, distance of BS 1 from 
the origin. A positive best response value indicates a location 
on the other side of the origin away from BS 1. Numerical 
computations indicate the existence of a unique symmetric 
equilibrium at —xi = X2 ~ 7.36. 

In Table U we compare the optimal location of the coopera- 
tive case and the equilibrium location of the non-cooperative 
case, as a function of a. We observe that at the non-cooperative 
equilibrium, the BSs are closer than at the cooperative op- 
timum, i.e., placements are more aggressive. An analytical 
proof of this observation remains an interesting open problem. 
In both cooperative and non-cooperative cases the distances 
decrease in tr and tend to a limit which is —xi = X2 = 5 
for the cooperative case and — xi ~ X2 ~ 4.06 for the non- 
cooperative case. 
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TABLE I 

Optimal cooperative and non-cooperative placements of BSs 

AS A function of a 



a 


0.1 0.4 1 2 40 


Optimum distance 
of BSs from 
(cooperative) 


8.658 7.745 6.435 5.591 5.002 


Equilibrium distance 
of BSs from 
(non-cooperative) 


8.10 6.95 5.50 4.667 4.09 



IV. CDMA; The Two-Frequencies Case 

A. SINR-equilibrium association 

We study the properties of the SINR-equilibrium partition 
and arrive at a numerical method to compute it. The case where 
BSs are collocated, i.e., xi ~ X2, is trivial. Let us assume that 
xi 7^ X2- Recall that the interference at BS j in the two- 
frequencies case is E{xj,Aj). As in Section UlI-AI define B 
to be the ath root of the ratio of the net interferences at the 
two BSs, i.e., 

n 1/" 

^ r-y I _ All I -t- 



E{xi,Ai)+a'^ 



Note that B e [B, 

-^m i n 



mm; -^max 



_E°{x2) + cr^ 
A location y E A2 if 

9{y~X2) 



E{x2,A2) + a^ 
, where 

1/q 

and Sniax = 



E^xi) + 0-2 



1/q 



Eix2,A2) 4 

^ (y- 



> 



9{v-xi) 



(9) 



X2f + l< {{y~x^f + l) B^ 

First consider B < 1. Proceeding exactly as in the proof of 
Proposition 13. II we get A2 to be the set of y such that a convex 
quadratic function of y is negative. Thus A2 is an interval and 
its complement Ai a union of at most two intervals. More 
precisely, the boundaries are given as follows. Define 

t{B) \xi - X2 



If t{B) < 1, A2 is empty, B^ = {E{xi,Ai) + a'^)/a'^ > 1, 
a contradiction. Thus t{B) > 1 and A2 is determined by the 
interval (as in the proof of Proposition I3.lt 

(51 (S),. 92(5)) 

X2 — xiB^ 



+ -v/r(i3)2-l,v/r(S)2-l 



This gives expressions for the end points of intervals that make 
up Ai and A2 in terms of B. In particular^ 



A2 



[-L,Lgi(S)J-L)U( [.92(5)1 L,i], 

[L9i(s)j-L,r92(i?)i 



L ■ 



Similar analysis could be done for i? > 1. In this case, it can 
be shown that 

Ai = [[g,{B)\^L,\92{B)]Ll 

A2 = [-L,lgi{B)\.L)Ui\g2{B)]L,L], 



*r™l-t = min{m, L}, 
min{max{m, —L}, L} 



max{m, —L}, and 



where gi{B), g2{B), and t{B) are the same as defined above. 
Finally for B — 1, ^ implies 



Xi + X2 



1 L 



-L 



and A-2 = 



Xi + X2 



To emphasize that Ai and A2 depend only on B, we write 
Ai{B) and A2{B). At SINR-equilibrium, therefore, B must 
be a solution to the fixed point equation 



B = 



£;(xi,Ai(S))+a2 



E{X2,A2{B)) 



11 a 



^■■F{B). 



(10) 



Theorem 4.1: The fixed point equation F{B) = B has a 
unique solution. 

Proof: We first prove that E{xi,Ai{B)) and 
E{x2, A2{B)) are continuous in B. By inspection, gi{B) 
and g2{B) are continuous for all B ^ 1. Straightforward 
calculations show that gi{B) — >■ [xi +X2)/2 and g2{B) 00 
as i? t 1, while gi{B) -00 and .92(5) {xi + X2)/2 
as B ^ 1. So the boundaries of Ai{B) and A2{B), after 
restriction to [— are continuous in B in [Bmin, ^max]- 
Thus E{xi,Ai{B)) and E{x2,A2{B)) (see (E)) and therefore 
F{B) are continuous functions of B in [Bmin , 5,„ax] ■ 

Next we show that F{B) is a decreasing function of B. Let 
B <1. We observe that if a ?/ satisfies ^ for some value of 
B, it will do so also for any larger value of B in [i?min, !]■ 
Thus, A2{B) is an increasing set function (order is specified 
by inclusion relation), and Ai{B) is a decreasing set function. 
For > 1, a similar argument shows that with B' = \/B, 
Ai{B') is increasing in B' , and A2{B') is decreasing in B' . So 
Ai{B) and A2{B) are decreasing and increasing set functions 
respectively, for B in [l,i3max] as well. Obviously E{-,A) 
strictly increases as A C [— L] increases. Hence F{B) is a 
decreasing function of B for B in [i?min, 5,„ax]- 

Finally, from ( fTol i. we see that F{B) e [i?min, 5i„ax]. Thus 
there is a unique fixed point of the equation F{B) ~ B. ■ 

As a simple example, consider the symmetric case when 
—xi = 2:2 7^ 0. It is easy to verify that B = \, and the 
unique equilibrium partition is ([— L, 0], (0, L\) if xi < and 
([0,L], [-L,0)) if xi >0. 

For an integer a > 1, ( fTOl i can be written as an implicit 
equation in B as 

B" = 

[arctanQ(L — a;i) — aictaj\a{\ g2{By\ l ^ xi) 

+ arctanQ,([gi(i3)J - l — xi) — arctanQ,(— L — I'l) + cr2j 

arctana([g2(5)]L - X2) - aicta.iia{[gi{B)\-L - X2) + 0-2 
We may numerically search for a B that solves the above 
equation through a suitably fine quantization of the specified 
interval. 

Remark 4.1: Suppose B* solves the fixed point equa- 
tion ( [Tot . B* <1 implies that, at equilibrium, BS 2 has more 
interference than BS 1, and thus A2 is a connected subset of 
[— L, L]. B* > \ implies that, at equilibrium, BS 1 has more 
interference than BS 2, and thus Ai is a connected subset of 
[-L,L]. 

Figure |8] illustrates the existence and uniqueness of fixed 
points. We set L = 10, the path loss exponent a = 2 and the 
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Fig. 9. Two frequencies case; SINR-equilibrium: Thu'esholds determining the 
cell boundaries as a function of the location of BS 2 for various locations of 
BS 1. 



Fig. 8. Two-frequencies case; F(B) vs B. Figure illustrates the existence and 
uniqueness of fixed points. 



noise parameter cr = 0.3. Two curves are shown: one corre- 
sponding to (xi, .T2) — (0, 10) and the other corresponding to 
{xi,X2) = (10,0). For both the cases F{B) is monotonically 
decreasing and cuts the Une y = x at unique points. In the 
first case, the fixed point is B ^ 1.393; at equiHbrium, BS 1 
has more interference than BS 2. In the diagonally opposite 
second case, the fixed point B = 0.718 (= 1/1.393). 

We also provide an algorithm to obtain the fixed point. 
Define d := \xi — X2I/2, /3,„in ;= V(P~fT — d and 
/3max V<PTT + d. Then, for B < I, t{B) > 1 implies 
B e (/3min, 1]- Similarly, for B > 1, t{B) > 1 implies that 
B e [1, /3max)- From an earlier discussion, the fixed point lies 
in (/3minj /3max)- It Can also be verified that F'{B) —00 
as B ^ /3min or /3max ^nd is finite otherwise. Thus a fixed 
point iteration of the form B^+i — F{Bn) (see (fTOl i) may 
not converge. Motivated by ifTSl . we propose a variant of this 
iteration which always converges to the desired fixed point. 
See Appendix ID] for the algorithm and its analysis. 

We illustrate the form of SINR-equilibrium via some numer- 
ical results. We set L = 10, the path loss exponent a = 2 and 
the noise parameter a = 0.3. We place BS 1 at one of the fixed 
locations xi where xi = —10, —5,— 2,0. For each of these, 
we vary the location of BS 2 from 0:2 = to 2:2 = 30 (see 
Figure |9]l. As in single frequency case, the equilibrium sets 
Aj have the form Ai = [6*1, 6*2], A2 = [-i,6'i) U {02, L] 
for xi = 0,-2, and Ai = [-L,02], A2 = (6*2, i] for 
Xi = —5,-10. The left (respectively right) plot depicts the 
threshold di (respectively 62) as a function of X2- 



at BS j. Then the sum utility satisfies the following: 



^ 1 

T-- 

^ 2 u. 



< 



< 



(Ul +U2)/2 



U2)/2- 



nax/2 + ^2 

E{LI2, [0,L] 



(11) 

(12) 
(13) 



E{L/2,[Q,L]) + a 

where ( fTTT i follows from Jensen's inequality because the 
function ii/(7i + cr^) is concave in u\ inequality (fTST i follows 
because the function 
with u. 



u/2 



—TTT-, — - IS monotone mcreasmg m u 
the maximum sum of received energies across 



any partition (not just SINR-equilibrium partitions). The last 
equality ( fT3] ) follows from Proposition IE. II in Appendix |E] 
The upper bound is independent of xi and X2, and is achieved 
when ~xi = X2 = L/2. The corresponding intervals indeed 
constitute an SINR-equiHbrium cell partition. ■ 
2) Two non-cooperating base stations: We now consider 
the hierarchical game where the BSs compete with each other 
keeping in mind their individual utilities as in Section IIII-B2I 



O.Sr 



0.4 



(Jl 0.3 
in 



■■ 0.2 




0.1 



B. Hierarchical equilibrium 

1 ) Two cooperating base stations: The goal here is to place 
the two BSs so that the sum utility is maximized. 

Proposition 4.1: The locations —xi = X2 = L/2 with 
SINR-equilibrium cell partition (^1,^2) = ( [-i, 0], (0,i] ) 
maximizes the sum utility. 

Proof: For a given pair of locations xi and X2, let 
(^1,^2) be the SINR-equilibrium cell partition. For conve- 
nience let Uj E{xj,Aj), j = 1,2, be the received power 



-% -5 5 10 15 20 

Fig. 10. Two frequencies case; non-cooperating BSs: Utility of BS 2 
as a function of its location when we position BS 1 at xi where xi = 
-2,-5,-8,-10. 

Figure [TO] has on the horizontal axis the location X2 of 
BS 2 and on the vertical axis the utility it achieves. The 
figure is obtained for L — 10, cr = 0.3, a = 2. There 
are four curves that correspond to four locations of BS 1: 
xi ^ —2, —5, —8, —10. From these curves, one can conclude 



9 



that the utility of BS 2 is quite robust to placement errors 
around the best response location, for the indicated values of 
BS 1 locations. Figure [TT| yields the best response for BS 2 
given BS I's placement. Given a BS 1 location, the higher 
interference cell and the equilibrium ratio B are first found 
as discussed in Section IIV-AI for each possible location of 
BS 2. Then the BS 2 location yielding the maximum utility 
is identified as the best response location and is plotted in the 
figure. 
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Fig. 11. Two frequencies case; non-cooperating BSs: The best response of 
BS 2 when BS 1 at a distance indicated by abscissa to the left of the origin. 
A positive best response indicates a location on the other side of BS 1. 

Numerical results indicate that there is a unique symmetric 
equiUbrium for the chosen parameters at — xi = X2 = 4.1. The 
corresponding SINR-equilibrium cell partition is Ai = [— L,0] 
and A2 ~ {0,L]. Note that any unilateral deviation will 
change the cell boundaries and yield lesser utility to the 
deviating BS. If the BSs were cooperative, the best locations 
are — xi = X2 = 5. However, this latter set of locations is 
not an equilibrium under competition. We also performed nu- 
merical computations for other values of noise variance, e.g., 
cr = 0.1, 1, 2. It is observed that the equilibrium BS locations 
are insensitive to the value of a. Yet again, we observe in our 
numerical examples that placements are more aggressive (i.e., 
BSs are closer to each other) under competition. 

V. SIC: The Single-Frequency Case 

We now extend our study to incorporate the effect of em- 
ploying SIC decoding by the BSs. For model description, see 
Appendix lAl In this section, we consider the single frequency 
case where mobiles connected to BS 1 cause interference at 
BS 2 and vice-versa. 



A. SINR-equilibrium association 

We first study the properties of the SINR-equilibrium parti- 
tion. The SINR density seen at BS j, for a mobile at y E Aj, 



can range from 



W17J)+^ '° E°{xj)-E{x-,Aj)+a''' 



, depending 



on the BS's decoding order We assume that each mobile 
first associates with a BS. The BSs then choose an arbitrary 
decoding order In the absence of a clear policy for choosing 



the decoding order at the BSs, we assume that the mobiles at 
y associate to BS 2 only if 



9{y~X2) 



> 



We may interpret this as an association where a mobile 
optimistically believes that it will be decoded last (at the BS 
with which it is associated) and therefore expects to see an 
SINR density (^E^ix ■)-%~{x'\ )+ai) Without loss 

of generaUty, relabef indices so that 

'E{xi,A2) + cr^ 



B 



l/c 



_E{x2,Ai) + a^ _ 
The above condition can be rewritten as 



< 1. 



{y-X2f + l< {{y-x^f + l)B\ 

where B"^ <1. 

The condition governing the structure of the SINR- 
equilibrium partition has the same form as the one in 
Section IIII-AI As before, the cell partition (^1,^2) = 
([-i,0], (0,L]) is an SINR-equilibrium with B = I m the 
symmetric case when —xi ~ X2- 

Numerical computations show that there can be more 
than one SINR-equilibrium partitions for a pair of BS lo- 
cations {xi^X2)- To illustrate this phenomenon, assume that 
all the mobiles are associated with BS 1. A mobile's antici- 
pated SINR densities at the two BSs would be and 
-E°(a:7)+o-^ respectively. It does not switch to BS 2 if 

g{y-x2) ^giy-xi) 



< 



E°{x2) + a^ 
Thus all mobiles keep their associations if 



g{y - X2) 



< 1 



E°ix2) 



yye[-L,L] 



giy-xi) 

This makes ([— L, L],0) an SINR-equilibrium partition. Sim- 



ilarly one can argue that 
equilibrium partition if 



g{y - xi) 



< 1 



(0, is also an SINR- 



giy-x2) 

If the noise variance is sufficiently small, both conditions 
can hold and hence both ([— L,i],0) and (0, [— L,i]) are 
SINR-equilibrium partitions. 

Remark 5.1: The above discussion illustrates an interesting 
capture phenomenon, which is of interest if the two BSs are 
placed sequentially, i.e., a Stackelberg game is played. The 
BS being placed first can judiciously place itself, to capture a 
majority of the mobiles. 

B. Hierarchical equilibrium 

BS 1 and BS 2 choose their respective locations coopera- 
tively. Each BS then employs SIC decoding for all mobiles 
in its cell. From the discussion in Appendix lAl the utility of 

independent of the decoding 



BSlisilog(l + ^f^^ 
order A similar expression is obtained for the utility of BS 2. 
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Two cooperating base stations: In the cooperative case, if 
cr^ 0, it is nearly optimal if all mobiles can associate to 
one BS. This is because if there is a non-zero population 
of mobiles connected to one BS, it generates a non-zero 
interference to the other BS. On the other hand, with all the 
mobiles associated to one of the BSs, say BS 1, the sum utility 



ilog( 



1 



oo when <t 



0. So BS 2 should be 



placed very far away so that its cell is nearly empty. Symmetric 
placements are therefore not optimal in general. However, as 
(7 — > oo, interference from the other BS no longer plays a 
role. As in the single mobile decoding scenario, BSs should 
be placed to maximize the total power collected from their 
respective cells. Thus symmetric locations — xi = •'^2 = ■§■ 
become optimal. 

Recall that there can be multiple SINR-equilibria for a given 
pair of BSs locations. Thus for the case of non-cooperating 
BSs, the competitive equilibria are not well defined. We do 
not pursue their study in this work. 



C. Pessimistically behaving mobiles 

In the absence of a clear policy for choosing the decoding 
order at the BSs, mobiles might also think pessimistically, 
i.e., each mobile makes an association decision assuming that 
it will be decoded first (at the BS with which it associates) 
and therefore expects to see an SINR density with 
BS j. This is the same SINR-density as observed in the case 
of single mobile decoding. Hence for given BS placements, 
SINR-equilibrium partition will be identical to that in the case 
of single mobile decoding (Section [III- Al l. 

For cooperating BSs, hierarchical equilibrium will be asym- 
metric in general as in the case of optimistic mobiles. When 
BSs are selfish, BS 1 is interested in optimizing 



1 



E{XI,A2)+(T^ 



■ log 



g°(a:i)+a-^ 
_ E{xi,Ai) 



g(ai,Ai) 



(14) 



which is the 



This is equivalent to optimizing E''{xi)+a^ 
utility of BS 1 if single mobile decoding is employed. In 
the single mobile decoding case numerical examples indicate 
that there is a hierarchical equilibrium with symmetric BS 
placements, say ~xi = xi = x*, and cell partition {Ai, A2) = 
{[-L,0], (0,L]) (Section liri-B2l ). Clearly, the same is a hier- 
archical equilibrium in the case of SIC decoding as well. For a 
BS, the aggregate equilibrium utilities corresponding to single 
mobile decoding and SIC decoding, are related as in (fT4l i. 

Remark 5.2: Having discussed the two extreme decoding 
order beliefs, we are naturally led to the following interesting 
problem. What decoding policy should a BS advertise in order 
to maximize its utility? Recall that given an association profile, 
a BS's utility does not depend on what decoding order it 
actually follows. But the advertisement will affect the SINR- 
equilibrium and thus the utility. We leave this as an open 
problem for future research. 



VI. SIC: The Two-Frequencies Case 

We now proceed to the two-frequencies case where the BSs 
employ SIC decoding. We give a complete characterization of 
both cooperative and competitive equilibria. 

A. SINR-equilibrium association 

In the two frequencies case, the SINR density seen at BS 7, 
for a mobile at j/ e Aj, can range from e(x^'a^)+(j'^ 
gCy-^-cj) ^ depending on the BS's decoding order As before, 
we assume that each mobile first associates with a BS. The 
BSs then choose an arbitrary decoding order In the absence 
of a clear policy for choosing the decoding order at the BSs, 
an optimistic mobile believes that it will be decoded last and 
therefore expects to see an SINR density of g{y — Xj)/a'^ with 
BS j. This being monotonically decreasing in the distance 
\y — Xj\, the mobile simply associates to the nearest BS. If the 
BSs are collocated, then either BS is chosen arbitrarily. Define 
V = (xi +X2)/2. Then, the equilibrium cell partition is Ai = 
[-L,v], A2 = if .Ti <X2, Ai = [v,L], A2 = [-L,v) 

if xi > X2, and an arbitrary choice at every y if xi = X2- 

B. Hierarchical equilibrium 

In the two frequencies case, BS j's utility with SIC decoding 
is i log(l + -^^^T^), j = 1, 2 (see AppendixlHi. 

1) Two cooperating base stations: In this case, the two BSs 
cooperate to maximize sum utility. 

Theorem 6.1: Consider the two-frequencies case with two 
cooperating BSs that employ SIC decoding. The BS locations 
that maximize sum throughput are —xi = X2 — L/2. 

Proof: Recall the notation used in the proof of Proposi- 
tion |4T| where Uj ~ E{xj,Aj). The sum throughput may be 
upper bounded as 
2 



Ui + U2 

2a2 



< log [I ^ 
= log(l + 



E{L/2,{0,L] 



where the first inequality follows from Jensen's inequality, 
while the second follows as in the proof of Proposition 14.11 
and the third follows from Proposition IE. II Finally, the upper 
bound is attained at — xi = .T2 = L/2 with cell partition 
( [— i,0], (0, L] ). This completes the proof. ■ 
Consider now a case where the BSs are constrained to be 
collocated at x. Recall that mobiles pick one or the other BS 
with equal probability, so that the power collected at each 

This 



BS is E°{x)/2 yielding a sum utiHty log [l + 
attains its maximum when E° does, which is at x = (see 
Proposition 12. 11 1. 

2) Two non-cooperating base stations: In this case, the two 
BSs play a non-cooperative game to maximize their respective 
utilities. Define a := 2^/". For a G [1, 00), we have a G (1, 4]. 
Recall that if the two BSs are not collocated, the cell boundary 
is (xi+X2)/2. Let rj{xi,X2) be the power collected by BS j. 
The utility of each BS is a monotone function of the power 
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collected, and we may therefore assume that BS j's goal is to 
maximize rj (x'l, X2). Let BR2(a;i) be the set of best responses 
of BS 2 to xi, i.e., 

BR2(.Ti) = argmaxr2(xi,a;2) 

Define BRi(x2) analogously. 

Lemma 6.1: (i) For all xi, BR2(.ti) C [— L,i]. An analo- 
gous conclusion holds for BRi(j;2). 

(ii) If (xi, a;2) is an equilibrium strategy profile, then xi,X2 G 

Proof: We consider the following three cases 

1) If Xi e (-00, -L), then r2{xi,-L) = E°{L). On the 
other hand if \x2\ > L, then r2{xi,X2) < E°{x2) < 
E°{L). Thus if an X2 ^ then X2 ^ BR2(xi). 

2) If xi e then 

r2(.Ti,-L) > r2(a;i,a;2), V.T2 e (-00, -i), 
andr2(a;i,L) > r2(a;i, 3:2), V.T2 G (i, 00). 

Thus if an X2 G (— 00, — i) U (L,oo), then X2 ^ 
BR2(a;i). 

3) If xi e (L, 00), then r2(xi,i) = E°{L). On the other 
hand if |x2| > L, then r2(xi,X2) < E°{x2) < E°{L). 
Thus if an X2 ^ [—£,£], then X2 ^ BR2(xi). 

Similar arguments hold for BRi(x2) also. 

The second part follows immediately once we recognize 
that (xi,X2) is an equilibrium strategy profile if and only if 
Xi G BRi(x2) and X2 G BR2(xi). ■ 

Thus, for equilibrium analysis, we only need to focus on 
a;i,X2 G [—L,L\. For xi,.T2 G [—L,L\, r2(xi,X2) is as 
given in Table HJ with a similar table for ri(xi,X2) of BS 1. 
Interestingly, the function r2(xi, ) as a function of X2 is 
discontinuous at X2 = Xi unless xi = 0. A similar observation 
holds for ri(-,X2). Next we show that in an equilibrium, the 

TABLE II 

Power received at BS 2: g [-L, L] 



X2 e 


r2(xi,X2) 


[-L,xi) 


arctana 2 } S'i'ctana {L + X2) 




E°{xi)/2 


{xi,L] 


arctariQ (L ~ X2) + arctaiia ^ ^2 -J^i ^ 



two BSs are placed on the opposite sides of the origin. 

Lemma 6.2: If xi < 0, then BR2(xi) > xi. Similarly, if 
xi > 0, then BR2(xi) < xi. Analogous conclusions hold for 
BRi(x2) also. 

Proof: Let xi < 0. If .xi < — L, the result follows from 
Lemma |6T| For — L < xi < 

r2(xi,X]^) > r2(xi,xi), 
and r2(xi, -X2) > r2(xi, X2), Vx2 < xi. 

See Table |ll] to verify the first claim above. The second one 
follows since mobiles associate to the nearest BS. These two 
imply that if X2 < xi, then .X2 ^ BR2(xi). This is the desired 
result. Similarly one can argue for xi > 0, and subsequently 
for BRi(x2). ■ 



Corollary 6.1: If (xi, X2) is an equilibrium strategy profile, 
then either (i) xi G [— L, 0],X2 G [0,L], or (ii) xi G 
[0,i],X2 G [-L,0]. 

Proof: If xi = 0, there is nothing to prove. First let 
xi G [— 0). The above lemma implies that X2 ^ xi, because 
X2 = BR2(xi) > x\. Moreover, .X2 <^ (.xi,0), because xi = 
BRi(x2) > X2 for X2 < 0. Thus X2 G [0,L]. Case (ii) is 
similarly handled. ■ 

We now characterize all the equilibria. 

Theorem 6.2: (i) For L < \/a — 1, there exists a unique 
equilibrium at xi = X2 = 0. 

(ii) For L > ^/a — 1, there exists a unique equilibrium (up to a 
permutation) at — xi = X2 = ( ^ ■ 
Proof: See Appendix |F] ■ 
Remark 6.1: (i) The equilibria locations do not depend 
on a. A similar insensitivity observation was made for the 
equilibria in single mobile decoding case (Section [IV-B2| |. 
(ii) Note that 



-L+^aL^ - (a - 1)2 



< 



L 



L 

T<2- 



a — 1 a — 1 y 

Again, as already seen in Sections IIII-B2I and IIV-B2I for the 
case of single mobile decoding, the competitive equilibrium 
locations of BSs are closer to each other than the optimal 
locations under cooperation. 



a = 1 






'<1=-1 

= -10 







-10 



10 



Fig. 12. SIC two frequencies case: r2(x\,X2) as a function of location X2 
of BS2 for thi'ee fixed locations of BSl: x'l = —1, —5, —10. Here L = 10 
and = 1 (left) and a = 2 (right). 

Figure [12] plots the received power at BS 2 given by 
r2(xi, X2) as a function of location X2 for three fixed locations 
ofBS 1: xi = -1,-5,-10. Here L = 10 and a = 1 (left) and 
a = 2 (right). Both satisfy L > — 1. A clear best response 
location X2 is seen for each fixed xi. As BS 1 approaches the 
origin, the best response location of BS 2 moves further away 
from the origin. 

3) Convergence to equilibrium: We consider the best re- 
sponse dynamics in which the location of each of the two 
BSs is sequentially adjusted. 

Theorem 6.3: Let L > ^/a — 1. Assume that BSs follow the 
best response dynamics to adjust their positions. Then, starting 
from arbitrary initial positions xj and .X2, the best response 
sequence converges to the unique equilibrium. 

Proof: See Appendix iGl ■ 

As an example, consider L = 10 and a = 2, i.e., a = 2. 
The equilibrium locations of the BSs are at a distance of 4.107 
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from the origin. Figure [T3] illustrates the fast convergence of 
dynamics from the starting locations ±5 to the equilibrium 
locations. 
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Fig. 13. SIC two frequencies case: non-cooperative BSs; Convergence of 
the best response dynamics for BS 2. 

Remark 6.2: When mobiles behave pessimistically, i.e., 
each mobile makes its association decision assuming that it 
will be decoded first (at the BS with which it associates), 
SINR-equilibrium and hierarchical equilibrium results are 
identical to those for the single mobile decoding in the two- 
frequencies case given in Section |IV] (See the discussion in 
Section IV-CI for a justification in the single-frequency case). 

VII. Conclusions 

We studied combined BS placements and mobile associa- 
tions in a game-setting where the utilities were determined 
by SINR criteria. We saw that the SINR-equilibrium cells 
exhibited non-monotonicity and non-convexity properties that 
are not seen in the classical location game problems. These 
unusual properties arise because the SINR density that de- 
termines association is a function of the distance between 
a mobile and the BS it is associated with and also the BS 
location. We studied hierarchical equilibria in the CDMA 
single-frequency and two-frequencies cases. We saw evidence 
(via simulations in the CDMA single-frequency case and via 
analysis in the other case) of a unique optimal pair of locations 
in the cooperative scenario. We also saw evidence of a unique 
equilibrium pair of locations (up to permutation) in the com- 
petitive scenario. For the SIC single-frequency case, we made 
some interesting observations. For the SIC two-frequencies 
case, we completely characterized the optimal cooperative 
locations and all competitive equilibria. Interestingly, in all 
scenarios considered, the BS locations are closer to each 
other in the competitive case than in the cooperative case, 
an observation whose proof has eluded us. 

Appendix A 
Propagation, path loss, and fluid models 

Propagation model: A mobile transmitter is modeled 
as a point source that radiates in two-dimensional space or 
three-dimensional space. The wavefronts emanating from the 



point source are circular (respectively, spherical in three- 
dimensional space). We assume that the far field model holds 
and that antenna couplings between neighboring transmitters 
and between transmitters and receiver are negligible, even in 
the limit as mobiles get closer to each other 

Path loss model: Under the far-field model for propa- 
gation in two dimensions with circular wavefronts, a receiver 
at a distance r from the point source and having aperture arc 
width s <C r will capture only s / (27rr) of the total transmitted 
power, so that propagation loss is proportional to 1/r. If there 
is further dissipation in the medium (analogous to shadowing 
and scattering of electromagnetic waves in three dimensions) 
we model the propagation loss as proportional to 1/r", where 
a > 1. The path loss model l/r" for three dimensional 
propagation with a > 2 is of course the standard one. 

Fluid model: Consider n mobiles located on a line at 
positions — L+jAy + -^, j = 0, 1, • • • , n— 1 with separation 
spacing Ay = — . We use the letter y to represent the discrete 
location for finite n, and the continuum location y E {—L, L) 
when n — )- 00. Each mobile has power Ap{y) = Ay, so 
that we may think of transmitted power density per unit 
distance dp/dy as 1 power unit per unit distance, and the 
total transmitted power as 2L power units. Consider the BS 
located at a; at a height of 1 unit from the line. The path loss 
for a mobile at y is g{y — x) = [l + (y — a;)^] (see ([T]i). 
The total received power at the BS, if all of these are in the 
same frequency band, is 



En{x) 



L-Ay/2 

E 

y=-L+Ay/2 
fL 

g{y~x)dy 

' -L 



g{y~x)Ay 



E"{x) 



where the limit is taken as n — s- 00. Similarly, the total received 
power from mobiles in a set ^4 C [— L, i] is 

En{x,A) = ^g{y-x)Ay 

yeA 

g{y~x)dy ^ E{x,A). 



E{x,A) was defined in ^ and E°{x) was defined as 
E{x, [—L,L]) immediately after. 

SINR Density and throughput: Let / C [—L,L] denote 
the set of locations that may be considered as interferer 
locations. The SINR is then 

CTMD t T\ 9{y~x)Ay 

SINR„(y,a;,/) = r-r- 

^ +2^yei9{v~ x)Ay 

As n ^ 00, the denominator tends to E{x,I) + a^, the 
numerator goes to 0, and the ratio 

SINR„(y,x,/) _^ g{y-x) 



Ay ' E{xJ) + a^' 

so that the latter may be thought of as SINR density (SINR per 
unit distance). Using Shannon's capacity formula for Gaussian 
channels, the data rate for a mobile at location y is 

ilog(l + SINR„(y,x,/)) « isiNR„(y,x-,/) 
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where the natural logarithm is employed and the unit of 
information is nats. (1 nat = l/(log2) bits w 1.44 bits). The 
aggregate throughput of mobiles in a set A C [— L,i] is 

V- 1. 



^2 

y&A 



SINR„(y,a;,/) 



1 



g{y - x) 



1 E{x,A) 



2 J A E{x, I) + cr2 2 E{x, I)+(J^' 

which is taken as the utility of a BS in the continuum case. 

SIC decoding: Let the interval A C [—L, L] denote a 
set of locations associated with the BS at x. Suppose that the 
BS employs SIC. An arbitrary decoding order is chosen and 
communicated with the transmitters. For concreteness, let us 
assume that mobiles are decoded in the decreasing order of y 
in A. Then all mobiles in A that are to the left of a given user 
at y will become interferers to y. The throughput for user at 
y E A is therefore 



g{y-x)Ay 




y'<y.y'e 



a9{v' - x)/^y 



y'<y,y'eA 



x)Ay 



1 



y'<y,y'eA 



g{y' - x)Ay 



Summing these up over discrete y G A, and passing to the 
limit, we get the aggregate throughput of all the mobiles in 
set A to be 

1, , J:yeA9iy-^)^y \ , 1, . E{x,Ay 

an expression that is used in Sections |V] and |VI] for the utility 
of BS. Note that this remains the sum utility regardless of the 
decoding order chosen at the BS. Of course, the data rates for 
each mobile will depend on its position in the decoding order. 
The sender and the receiver should agree on this data rate and 
employ an appropriate code. 

Discussion: It should be noted that the two-dimensional 
propagation (when a G [1,2)) and our treatment of mobiles 
as fluid particles on a line are merely caricatures of real 
life propagation models. The purpose of their study is to get 
a qualitative feel for what one might expect in the three- 
dimensional propagation model with mobiles distributed in a 
plane and receiver antennas placed at a height from the plane. 
See also the extension in Appendix |H1 

Appendix B 
Proof of Proposition I2.1I 

That E° is an even function, is obvious from (|6]l. To see 
the monotonicity, for x > 0, write 

E"{0) - E°{-x) = f g{y) dy - f ^ g{y) dy (15) 



L J-L+x 

— L-\-x pL-{-x 

giy) dy- g{y) dy 

-L JL 

[g{y - L) - g{y + L)] dy (16) 



where (flSl l follows from (|3), and (fT&t via a change of variable 
y — L ^ y in the first integral and y + L -s— y in the second. 
The integrand in (fT6] l is positive for y E [0,x]. This proves 
the monotonicity. 



Appendix C 

cdma single-frequency case: structure of 
cooperative hierarchical equilibria 

a) In an optimal configuration the two BSs are placed 
on the opposite sides of the origin: Assume that the two BSs 
are placed at xi,X2 such that xi^X2 < 0. The configuration 
xi = a;2 < is outperformed by xi — X2 = 0, and hence 
cannot be optimal. Next consider xi < X2 < Q without loss 
of generality. Under equihbrium, the cell of BS 2 is a 
convex set. Say A2 = [—a,h] C [— L, L]. The following are 
observations that are easily seen. 

• If a < 6, the configuration {xi, —X2, A2 = [—a, b],Ai = 
[—L, L]\[~a, b]} has better joint performance. 

> If a > 6, the configuration {xi, — 2:2, A2 ~ [—6, a],Ai — 
[— L, L]\[—b, a]} has better joint performance. 

• If a = 6, the configuration {xi, — X2, A2 = [—a, b],Ai = 
[—L,L]\[—a,b]} performs equally well. 

In each case, the SINR-equilibrium of the modification can 
only have better joint performance. This proves the claim. 

b) In an optimal configuration, the cells of the two BSs 
are convex: Without loss of generality, we may assume xi < 
< X2 and that |xi| > |x2| > 0. Also assume that A2 — 
[—a, 6] such that a,b < L so that Ai = [—L, — a)U(6, L] is not 
convex (see Figure O. But then the configuration {xi,.T2 + 
L-b,Ai^ [-L,-a+{L~b)),A2 = [-a + {L-b),L]} (see 
Figure [T4] | strictly outperforms the assumed one. 



-L xi -a + {L-b) X2 + (L-b) L 

Fig. 14. In the above two figures A2 = [—a, b] and [—a + (L — b), L], re- 
spectively. The top one can not be an optimal configuration; it is outperformed 
by the bottom one. 



Indeed, BS 2's utility remains the same, and a few of 
BS 1 mobiles have moved closer, keeping the same associated 
population size. This yields a contradiction and proves the 
convexity claim. 

The above observation implies that in a hierarchical optimal 
configuration, the cells are of the form [— L,a] and (a, L]. 
Furthermore, the optimal configuration is also yielded by the 
following hierarchical problem. For a given partition, i.e., a 
fixed a, the BSs are placed to optimize throughputs of the 
respective cells. However, the partition is made to maximize 
the sum utility of both the BSs while keeping in mind their 
subsequent selfish placements. 
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Appendix D 
The Relaxed Fixed Point Iteration 

The algorithm: (i) First we determine the starting point. 
We aim at finding a point Bq such that both Bq and F{Bo) 
lie in (/^mim /^max)- Recall that both B and F{B) belong to 
[-Bmin , ^max] and F{B) is decreasing in B. Also, we can 
easily verify that F{B) = S,nax for all B iifl [Bn^n, (3min\, 
and F{B) ~ -Bmin for all B in [/3max, ^max]- A suitable choice 
of Bq depends on the order of the quantities Bmin, Bmax, /3min 
and ;5inax- The following are the four possible orders. (Note 
that the fixed point lies in (Bmin, Bmax) n (/3min, /3max), so 
the two intervals must intersect.) 

1) If Anin < Binin < Bmax < Anax, choOSC a Bq E 

2) If Bmin < Pmin < Bmax < /?max, there are two 
possibilities. 

a) If i^(Bmax) > /3n 

b) If F(Bmax) 
{l^min^F (/5min))- 

3) If /3min ^ Bmin ^ /^max 

possibilities. 

a) If F(Bmin) < /^rnax, choOSC Bq G (Bmin, /3max)- 

b) If F(Bmin) > /3max, choOSC Bq G 

(/^max); /^max)- 

4) If Bmin < /3min < ;3max < Bmax, choOSC a Bq G 
(i^-H/3max), i^-l(/3min)). 

Set B = F{Bo). Our selection procedure ensures that Bo and 
B are in (/3min, Anax)- If B = Bq, it is the desired fixed 
point. So, assume B > Bq; the case B < Bq is handled 
similarly. From our earlier discussion F'{B) is bounded over 
[Bo,B] C (/3niin, /3max), 1-6., there is a B < oo such that 
|F'(B)| < B for all B €[Bo,B]. Now, we focus only on this 
interval. 

(ii) Choose a 7 < and define 

G(B) := ^F{B) + (1 - 7)B, B G [Bq, B]. (17) 

(iii) Now iterate as B„+i = G{Bn), n — 0,1, ... . 
Proposition D.l: The iterates converge to the fixed point of 

F{B) = B. 

Proof: We prove that iterates as defined in Step 3) form 
a nondecreasing sequence. From (fTTT i. it suffices to show that 
F{Bn) > B„ for all n. We show this inductively. Since B > 
Bo, the claim is true for n = 0. Assume -F(B„) > B„ for 
some n. From (fTTj i. we conclude that B„+i > B„. Since 
|F'(B)| < D, we get 



B(B„) - B(B„+i) < B(B„+i - B„ 



(18) 



Moreover, our choice of 7 (see Step (ii)) together with (fTTI i 
ensures that 

B(B„+i — B„) < F{B„) — B„+i. 

Combining the above two inequalities we get F{B„+i) > 
Bn+i- This completes the induction step. 

Note that ( fTSl l requires B„,B„+i < B which can also 
be shown inductively, as follows. Since Bi is a convex 

5[a,fel = if a > fe. 



combination of Bo and B := F{Bq), Bi < B. Now assume 
Bo < Bi < . . . B„ < B. This immediately implies that 

B„+i = 7F(B„) + (1 - 7)B„ 
< 7F(Bo) + (1 - 7)B 
= B. 

Since the sequence B„,7i = 0,1,... is bounded and non- 
decreasing, it converges. Finally, since G{B) is a continuous 
function of B, the limit point B* is a fixed point of G(B) = B. 
Substituting in ([TtIi, we get F{B*) = B*. Therefore B* is the 
desired fixed point. ■ 
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Fig. 15. Two-frequencies case; F(B) vs B. The four subplots illustrate the four 
possible scenarios as described in Step 1 of the relaxed fixed point iteration. 



The four subplots in Figure [15] illustrate the four possible 
scenarios described in Step 1 above. We set L = 10, a = 2 
and cr = 0.3. Each subplot shows the corresponding locations 
of BSs. Note that, in all the cases, d — \xi — X2\/'2 = 2.5 
which also fixes (/3min, /3max) — (0.1926,5.1926). Bmin and 
Bmax corresponding to the four subplots are as in Table |III] 



TABLE III 

^min AND Bmax IN THE FOUR SUBPLOTS OF FIGURe[T5] 





Xl =15 Xl = 10 Xl = 5 XI = 

X2 = 10 X2 = 5 X2 = 10 X2 = 5 


-^min 


0.2364 0.1741 0.2364 0.1741 




1.6580 4.2306 5.7423 5.8045 



Figure [16] shows convergence of the relaxed fixed point 
iteration via an example: B = 10,a = 2 and a — 0.3; BS 1 
and BS 2 locations are set to be —20 and —15 respectively. 
Thus d = 2.5, and (/3mi„,Anax) = (0.1926,5.1926). Also 
Bmin = 0.6031 while Bmax = 1.3180. We set ^ = 0.1 and 
start the iteration with Bq = 1.3. Iterations quickly converge 
to the unique fixed point 0.726. 
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Fig. 16. Two frequencies case; convergence of the relaxed fixed point 
iteration. 



Appendix E 
Two-frequencies case: sum received power 

Proposition E.l: Let xi < xi and let v — (x\ + X2)/'2, 
denote the mid-point. Then following results hold. 

(i) Let (^1,^2) denote a partition of [— L,i]. Then 

max [E{xi,Ai) + E[x2,A2)\ 
(.A1.A2) 

< E{xi,[-L,v]) + E{x2,{v,L\). 

(ii) Furthermore, E{xi, [--Ljv]) + E{x2, {v,L]) is maximized 
at ~xi = X2 ~ L/2. The cell partition in this case is [— i, 0] 
and (0,L]. 

Proof: The first statement is obvious once we write out 
the integrals and recognize that the integrand is non-negative, 
symmetric, and g{y) is decreasing in For the same reason, 
we may upper bound the sum in the second statement, by 
E{xi,Ii) + E{x2, 12) where Ii is an interval of length L + v 
centered at xi and I2 is an interval of length L — w centered at 
2^2 ■ {{Iijli) may not be a partition of [— L, L]). Constraining 
the sum of interval lengths to be 2L, the upper bound is 
further maximized when the intervals are of equal length L. 
But this upper bound is achieved when —Xi = X2 = L/2. 
The corresponding intervals [— L, 0] and (0, L] constitute a 
cell partition, as required. This concludes the proof. ■ 

Appendix F 
Proof of Theorem I6.2I 

Let {xi,X2) be an equilibrium. We may assume, without 
loss of generality, that xi E [— L,0],.T2 G [0,L] (see Corol- 
larv l6.lt . Differentiating r2{xi,X2) with respect to X2, we get 

dr2{xi,X2) 



dx2 



2^ 



X2 - Xi 



g{L-x2) (19) 



The rest of the proof is divided into two parts. 

Part 1: We first check whether xi = 2:2 = is an 
equilibrium. Note that r2(0,X2) is an even, continuous and 
differentiable function of X2- From (fT9] l 

9r2 (0,0:2) 



dx2 







2Lx2 + - (a - 1) = 0, 



which implies 



:±v/^ + (a-l)(l-f) 



X2 



2L ■ 



, ae(l,4) 
a = 4 

Now, it is easy to check that ^''"g^";''"^ < for all X2 e (0, L], 
if L < Va^~T"- Thus, BR2(0) = 0. Analogously, BRi(O) = 
0. So, (0,0) is an equilibrium, provided L < ^/a — 1. For 

1 



L > y/a- 1, 

dr2{0,X2) 



dx2 



2:2=0+ 



2 
1 

^ 2 
> 0. 



m-9iL) 

1 



(l+L2)a/2 



Hence BR2(0) ^ 0. Thus (0, 0) can not be an equilibrium. 

Part 2: Next, let {xi,X2) 7^ (0, 0). Then, by assumption, 
we have X2 > xi. From ( fT9] l, a best response X2 to BS I's 
location xi should satisfy 



(L — X2)'^ + 1 = a 



1 



X2 — Xl 



(20) 



Similarly, a best response xi to BS 2's location X2 should 
satisfy 



(L + xi)^ + l = a 



1 



X2 — Xl 



(21) 



Combining ( l20l i and (l2Tl i we get 

{Xi + X2){X2 - Xl - 2L) = 

If X2 — Xl = 2L, the only feasible candidate is (a;i,a;2) = 
{—L,L). But, this contradicts (I2TI 1. Hence X2 = —xi. Again, 
from (I2TI 1. we have 



(L + Xif + 1^ a{l + xl) 



which implies 



Xl 



L±^aL^ - (a - 1)2 
a- 1 



Since xi G [— i,0], we must have aL^ ^ {(^ ^ 1)^ > 
so there will be no feasible solution for L < ^/a — 1. 
Furthermore, 



Xl 



-X2 



i-L + yZc^^T^^IF) 
a — i V / 



is the unique feasible candidate for L > \/a — 1. Finally, let 

X* := ^ \JO'L'^ - {a - 1)^)- We now show that 



Xl = —X2 = X* is indeed an equilibrium if i > \/a — I 
To see this, note that 

d'^r2{x*,X2) 



dxl 



1 



^-g'{-x*)+g'{L + x*)<Q, 

xi = — x' 

because x* < and L + .t* > 0. This assures that — x* is 
indeed a maximum point, i.e., —x* = BR2(x*). Similarly, one 
can verify that x* — BRi(— x*). 

Parts 1 and 2 together complete the proof. 
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Appendix G 
Proof of Theorem |6.3I 
Lemma 16.11 implies that it is sufficient to show the result 
for the case when, x°,X2 G [— i,L]. Assume xi < 0. 
Solving ( l20l i. we get 



X2^<, ^^^^ , 4-a ' (22) 



I. 2 2(L-xi)' u, — 1 

For a € (1,4), the solution to (|20] | with the positive sign 
leads to X2 > L. We therefore discard it, because BR2(xi) G 
[0,i]. The other candidate is the X2 in (|22] | which satisfies 
xi < X2 < L. This further implies that 

d'^r2{xi,X2) 1 



79 



X2 - Xi 



+ .g'(L-a;2) < 0. 



Hence X2 as given in ( |22] | is indeed BR2(.ti). Now, it can 
also be seen that < ^^g^^"^'^ < 1 - e where e > depends 
on a. An analogous result holds for xi > also. Analogous 

results also hold for ^M^ii£2) Thus a small change in the 

01:2 

position of a BS causes an even smaller change (in the same 
direction) in the position of the other BS. The best responses 
thus constitute a contraction map and the dynamics converges 
to the equihbrium. 

Appendix H 

cdma single-frequency case: extension to two 
Dimensions 

A large number of mobiles are located uniformly over the 
two dimensional plain. Two BSs are placed at (a;i,yi,l) 
and (.T2, 2/2,1) respectively. The BSs operate on the same 
frequency band. The path loss model is as before. The SINR- 
equilibrium association can be defined in a similar way as in 
Definition 12.11 Here, we provide closed form expressions for 
cell boundaries in the equilibrium. As in Section Illl-AI we 
define 

^ ^f E°{x,,v,) + a ^^"'^ 
" \E°{x2.y2) + o 
When the two BSs are collocated, every association profile is 
an SlNR-equilibrium. If (xi,yi) 7^ (3:2,2/2) and still i?^ = 1, 
there is unique SlNR-equilibrium association in which mobiles 
associate with the BS which is closer To study the asymmetric 
scenarios, without loss of generality assume that < 1, i.e., 
the interference at BS 2, E°{x2, 2/2), is more than that at BS 1, 
E°{xi,yi). Mobiles at (x, y) connect to BS 2 if they do not 
have a lower SINR density at (0:2, j/2), i.e., 

[(a;-.T2)" + (2/-;/2)' + l]~°/' 
E''{x2,y2) + (7^ 
> [(x-a;i)^ + (y-yi)^ + l]-°/^ 
E°{xi,yi) + (7^ 

which is equivalent to 

(x - X2f + {y- 2/2)' + 1 < ((x - xif + iy~yif + l) bI 

This inequality straightforwardly yields that the set connecting 
to BS 2 is nonempty only if 



Furthermore, in case of a strict inequality above, the set is a 
disc with center 

X2 - xiBl 2/2 - yiBl 



l-Bl 



l-Bl 



and radius Vt^ — 1. 

Note that in the one-dimensional case, the BS with higher 
interference had an interval as its cell. The complement of this 
interval joined the second BS. The analog of the interval-type 
cell in two dimensions is a disc-type cell. Further, observe 
that we did not make any assumptions on the population 
density. The cell partition is always a disc and its complement. 
However, the population density does affect the interference 
and therefore the actual parameters of the cell partition. 

References 

[1] F. Plastria, "Static competitive facility location: An overview of opti- 
mization approaches," European Journal of Operational Research, vol. 
129, pp. 461^70, March 2001. 

[2] J. J. Gabszewicz and J. F. Thisse, Location, ser. Handbook of Game 
Theory with Economic Apphcations, R. Aumann and S. Halt, Eds. 
North Holland: Elsevier, 2003, vol. 1. 

[3] V. Mazalov and M. Sakaguchi, "Location game on the plane," Interna- 
tional Game Theory Review, 2003. 

[4] H. Hotelling, "Stability in competition," Economic Journal, vol. 39, pp. 
41-57, 1929. 

[5] A. Aram, C. Singh, S. Sarkar, and A. Kumai; "Cooperative profit 
sharing in coalition based resource allocation in wireless networks," in 
Proceedings of IEEE Infocom 2009, Rio de Janeiro, Brazil, April 2009. 

[6] A. Aram, S. Sarkar, and C. Singh, "Provider customer coalition games," 
in Proceedings of Gamenets 2009, Istanbul, Turkey, May 2009. 

[7] E. Altman, A. Kumar, C. Singh, and R. Sundaresan, "Spatial sinr 
games combining base station placement and mobile association," in 
Proceedings of IEEE Infocom 2009, Rio de Janeiro, Brazil, April 2009. 

[8] S. Ramanath, E. Altman, V. Kumar, V. Kavitha, and L. Thomas, "Fair 
assignment of base stations in cellular networks," in 22nd World Wireless 
Research Eorum (WWRF) Conference, Paris, France, May 2009. 

[9] A. Silva, E. Altman. M. Debbah, H. Tembine, and C. Jimene, "Optimal 
mobile association on hybrid networks: centralized and decentralized 
case," arXiv:0911.0257vl. 
[10] G. Kasbekar, E. Altman, and S. Sarkar, "A hierarchical spatial game 
over licensed resources," in Proceedings of Gamenets 2009, Istanbul, 
Turkey, May 2009. 

[11] J. G. Wardrop, "Some theoretical aspects of road traffic research," in 
Proceedings of the Institute of Civil Engineers, vol. 1, no. 36, 1952, pp. 
325-378. 

[12] W. H. Sandholm, "Potential games with continuous player sets," Games 
and Economic Behavior, vol. 97, no. 1, pp. 81-108, March 2001. 

[13] S. C. Dafermos, "The traffic assignment problem for multiclass-user 
transportation networks," Transportation Science, vol. 6, no. 1, pp. 73- 
87, Febmary 1972. 

[14] C. I. Chek and J. Jose B. Cruz, "Stackelberg solution for two-person 
games with biased information patterns," IEEE Transactions on Auto- 
matic Control, vol. AC-17, no. 6, pp. 791-798, December 1972. 

[15] A. Kumar, E. Altman, D. Miorandi, and M. Goyal, "New insights from a 
fixed point analysis of single cell ieee 802.1 1 wlans," IEEE Transactions 
on Networking, vol. 15, no. 3, pp. 588-601. June 2007. 



V {Xl - X2f + (2/1 - V2i^Ba_ 

l-Bl 



> 1. 



5.5 



5<^ 



CO 
OQ 

H — 

o 

c 4.5 



OS 

o 
o 



o o O O O <) 



3.5 







4 6 
Iteration number 



8 



10 



